On the excluded space in applications of Feshbach projection 

formalism 

S. Karatae; lidifl 

Department of Physics and Electronics, Rhodes University, 
P.O. Box 94, Grahamstown, 6I4O, South Africa. 



K. AmofiO 

School of Physics, University of Melbourne, Victoria 3010, Australia 

(Dated: February 2, 2008) 

Abstract 

Various model applications in nuclear structure and reactions have been formulated starting with 
the Feshbach projection formalism. In recent studies a truncated excluded space has been enu- 
merated to facilitate calculation and identify a convergence in expansions within that truncation. 
However, the effect of any remainder must be addressed before results from such can be considered 
physical. 
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Interest in the Feshbach projection formalism as a way to unify reaction and structure 
theory has increased in recent years with studies such as those of the no-core shell model 
(SM), like those made using G- matrix interactions py], and of the continuum shell model j2[. 
With reaction theory, the formalism lies at the heart of the recent multi-channel algebraic 
scattering (MCAS) theory [3] though applications of that method j3| so far assume that all 
important channels have been included in the V space. However, the formalism has long 
been used to make the framework of the optical model for elastic scattering jl, [6|. In that 
use, all channels other than those of the elastic scattering are taken to form the excluded 
space Q while the elastic channels span the space supposedly treated exactly. The dual space 
formalism is elegant but it is almost totally impractical due to the (usually large) number 
of reaction channels to be considered in defining terms such as the dynamic polarization 
potential. Often then, either a total phenomenological approach or a restricted selection of 
terms of importance in the excluded space need be made to facilitate evaluations. While the 
foregoing centered on reaction theory, such is also the case when the formalism is applied to 
consider properties of an isolated nucleus as with the G-matrix SM of Navratil and Barrett [l| 
and with the Continuum Shell Model (CSM) of Volya and Zelelevinsky [3]. The question 
arises as to the convergence of solutions obtained in the truncated (Q) space to those of the 
actual complete space. 

The dual space form of the Feshbach theory [s], 0] fragments the Hilbert space {\&} into 
subspaces V and Q by the action of projection operators P and Q respectively on the space 
spanned by the eigenfunctions \E' of the full Hamiltonian. The subspace V is spanned by 
the functions while that of Q is spanned by . The assumption is that the solutions 
within P space can be evaluated while that of the full space V + Q cannot. As projection 
operators, P and Q satisfy the conditions 

p2 = p. Q^ = Q- PQ = QP = 
P + Q = l. (1) 

Then, with the notation Hxy = XHY, where X, Y are any combination of the operators 
P, Q, the complete Schrodinger equation can be recast as 

{E - Hpp) P^ = HpqQ^ 

{E-Hqq)Q^ = HqpP^. (2) 

Using the second of these to define Q"^ in terms of P\I' then defines the V space equation 
to be solved, namely 



E — Hpp — HpQ— — — zHqp 
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P^ = 0. (3) 



Thus the contribution to the full Hamiltonian from coupling to the excluded states (the 
so-called "doorway" states) is defined as the additional term 

AH = HpQ \ Hqp. (4) 

[E - Hqq) 

If one restricts consideration first to nucleon-nucleus (A^^) scattering, where the V space 
may be taken to define only elastic scattering, then AH is the dynamic polarization po- 
tential aspect of the optical potential. It is an energy-dependent problem to evaluate. At 
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low energies, as evidenced by MCAS studies Q, and for light masses particularly, channel 
coupling of the nucleon to excited states of the target are the essential elements. The DPP 
that results is highly non-local, complex, and energy- and angular momentum dependent. 
One consequence is that phenomenological local potentials which are often used are not 
physically justified. For such cases it seems that the Q space may be enumerated with cou- 
pling of the incident nucleon to relatively few states of the target nucleus. At higher incident 
energies, values for which the giant resonances of the target may be excited, past studies 
have noted their influence in scattering as doorway states in second order processes jl]. 
Other channels which reflect in those studies require residual complex (phenomenological) 
optical potentials; they should be nonlocal. When the available energy coincides with the 
target in its continuum, a successful way to enumerate Q space effects is to adapt the KMT 
theory That has been done, for example, with the so-called (^-folding optical potential 
(lo| . For those energies, one cannot enumerate Q space in any detail but its effects in many 
applications to date [lO|, lul seem encompassed in the medium and Pauli blocking effects 
defined in the infinite matter (7-matrices upon which the (^-folding model is based. 

A truncated Hilbert space is also the key feature of the shell model for nuclear structure; 
being essential for practical calculation. One may define the V space as that in which the 
SM interaction is defined. The Q space then encompasses all of the higher hu admixtures 
lying outside and which can have effects due to long range correlations. It is well known 
that by restricting the V space to be just the Ohio for valence nucleons requires polarization 
charges, typically of 0.5e, to match the associated model results to observation. Such is a 
reflection of higher huj correlations. However, by increasing the basis size, e.g. to encompass 
a complete (0 -|- 2)huj space, then there remain yet higher order correlations that may still 
play an important role. Whatever the space enlargement with a SM, one still has missing 
pieces by definition. 

The G-matrix SM model interaction of Navratil and Barrett [l| as used later 
considers coupling up to three-body correlations. In that no-core SM, the Q space is assumed 
to be made of states including only those 2- and 3-body correlations above those included 
in the V space. However, then P + Q 1 hj construction. Higher order correlations may 
become important as the mass of the nucleus increases, and with them, the probability of 
forming large clusters. ^^C is a case in point. First there is the success of the a-cluster model 
in explaining the super- deformed 0+ state at 7.6 MeV, which suggests 4-particle correlations 
at least must be included in the SM. Second, it was observed that the excitation energy 
of the 2^ state diverged from the experimental value as the model space was increased. 
Finally the predicted B{E2; 2^ 0^) is not in good enough agreement with the measured 
value. Even including a three nucleon force, as was done recently [3], did not make much 
improvement. By contrast, when using the (fitted) MK3W interaction in the (0 -|- 2)hu} 
model space, very good results for the spectrum and for the electron scattering form factors 
were found for ^^C [l5|. The action of the fitting in determining the interaction implicitly 
includes higher-order correlations not considered in the those ab initio interactions. While 
the fitted interactions are not ab initio, the use of them does illustrate the effect of the Q 
space. 

The GSM calculations [2,] allow coupling to the continuum via one- and two-particle 
excitations. They consider the effects of coupling to the continuum (their restricted Q 
space) above the standard shell model {V space). By so doing they ascribe widths to 
resonance states without changing the energy of that state from what was found from a 
V space calculation. However, coupling to the Q space necessarily requires that there be 
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contributions from such coupling also to the real part of energy eigenvalues. The effects of 
higher order correlations remains also unknown for this model. 

To consider such limitations made upon the actual Q space of any problem, we make a 
three space development of the Feshbach formalism. The Hilbert space is divided now into 
the V, a reduced, enumerable, Q, and a new remainder TZ spaces. There are three associated 
projection operators, P, Q, and R satisfying P + Q + R = 1. In this case the Schrodinger 
equation can be cast into the three projected, coupled equations 



{E - Hpp) = {HpqQm + HprR^) 
{E - Hqq) = {HQpPm + HqrR^!) 
{E - Hup) R^ = {HppP^ + //rqQ*) . 



(5) 



These equations reduce to those of Eq. ([2]) in the limit that couplings to TZ space can be 
neglected. Such seem to be the case for NA scattering at low and intermediate energies. 
Rearranging Eqs. gives 



P^ 



(E-H 



QQ) 



[HqpP^ + HqrR^] 
[HppP^ + HpQQ^il] 



{E - Hpp) 

Using the second equation for in the first, and rearranging, gives 
1 „ 1 



(6) 



{E-H 



-H 



QQ) 



QR 



{E-H^ 



-H. 



RR) 



RQ 



1 



Hqp + 



1 



QR 



{E-H, 



-H, 



RP 



RR) 



{E - Hqq 

Likewise one can eliminate Q"^ from the Eqs. ([6]) and rearrange the result to find 
1 „ 1 



P^. (7) 



n -1 



1 - 



(E - Hpp) 



-H. 



RQ 



(E - Hqq) 



-H 



QR 



1 



{E-H^ 



RR) 



Hpp + H 



RQ 



{E-H, 



-H 



QQ) 



QP 



P^. (8) 



Then substituting Eqs. ^ and (IE]) into the first of Eqs. (Ill) gives the V space equations 

1 „ 1 



E — Hpp — H 



iPQ 



n -1 



1 - 



(P - Hqq 
1 



-H 



QR 



{E-H^ 



-H. 



RR) 



RQ 



(E- 
Hpp 



H 



QQ) 



Hqp + Hqp 



-H 



1 



(P — Hpp] 
1 

(P — Hpp) 



(P - Hpp] 
1 



RQ (TP TT 

Hpp + HpQ 



Hpp 

^Hqp 
1 



{E-H, 



-H 



QQ) 



QP 



P^i = 0. (9) 
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Again for situations where coupling to IZ space can be ignored, [HpR = Hqr = 0), this 
reduces to the form given in Eq. ([2]). If expansion is restricted to third order in couphngs, 
each inverse bracket term in Eq. (Q becomes the unit operator so that a reduced V space 
equation is 

E - Hpp - HpQ— — tHqp - HpR— — -Hrp 



Thus the couphng to excluded space of a problem may be via coupling to either the Q 
or TZ subspaces at second order in the expansion as well as via the chains V —>■ — »• 
TZ{Q) —>■ V in third order. The states in TZ space more than play the role of hallway states in 
past studies, for we now allow for states therein that directly connect by couplings omitted 
in the enumerated Q space. To the extent that coupling to those states in TZ space is not 
negligible, such have effect even at second order in expansions. Thus, unless one can ensure 
all important couplings are collected within an enumerated Q space, then no matter how 
much, and to what order, the effects of the enumerated Q space are taken, errors arise at 
the second order level. Results of using the G-matrix interaction developed by Navratil and 
Barrett ll|L in large space no-core shell model calculations for systems (such as ^Li [IG] and 
^^C which exhibit features of a-clustering, and so of four-body correlations, give 

evidence of such omission. That approach does show convergence to the enumerated Q space 
complete values, but those values remain far from the empirical ones. The problem starts at 
second order. Indeed, neglect of the higher order correlations in that G-matrix interaction 
has already been identified as a potential problem [l6|. In contrast, the earlier interactions 



developed by Zheng et al. [17| seem to accommodate such higher order effects through the 
way they form the G-matrices, as the values of the B{E2) in ^^C change from 5.45 to 6.9 
e^fm^ for evaluations in a Ohw and a (complete) no-core (0 + 2)hLj shell model, respectively. 
The results, obtained using the OXBASH shell model code [l8|, are reasonable given that 



the measured value is 7.77 e^fm^. But we note that in other studies 19|, made using a 
standard shell model, problems of only including 2huj excitations were found. In particular, 
the neglect of excitations above the 2huj level forces a much deeper binding energy for the 
ground state. Addition of the Ahw excitations, in part, resolved the problem, as has been 
illustrated in the case of ^^O 20 . 



While a formalism with which the effects of an enumerable Q space can be accounted may 
be aesthetically pleasing, it must be remembered at what level of evaluation problems can 
occur when seeking to apply same to reality. Certainly the approach is built upon a good 
theoretical base and is an excellent mathematical exercise. But for use in considering real 
systems, it is not necessarily nor sufficiently complete even at second order. That may be 
particularly problematic for deformed nuclei. The object of structure and reaction theory 
is to determine as much information as possible regarding real nuclei and in approaches 
built upon a Feshbach formalism, all strong coupling features need to be accounted at 
whatever order they appear. In principle there may no solution to this problem built from 
theory made using first principles and one may only have the specification of effective forces 
(phenomenological or otherwise) to enable structure and reactions to be dealt with in any 
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form of unified way. It would seem that the most fruitful way to proceed is to seek the 
dominant states in the spaces to which to couple. 
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